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Outline 


• Motivation 

- Add reflected solar radiance simulation in the exist fast infrared 
PCRTM 

- Increase the computation speed of RS radiance simulation for cloudy 
atmospheres while keep high accuracy 

• LUTs for Multiple Scattering of Cloud and Aerosol 

- Generation 

• Stream Dependence in DISORT Simulation 

• Remove Stream Dependent Bad Points 

- Compression 

• From 3.3 GB to 0.024 GB 

- Rebuilt Accuracy 

□ ABRDF < 10 -4 

• TOA Reflectance Simulation 

- Accuracy 

□ ARTOA < 1 0- 3 

- Speed 

• Five orders faster than 52 -stream DISORT 

• PCRTM_SOLAR is faster than the exist infrared PCRTM 

• Example of Application on IASI Spectral Data 
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Motivation 



• RTM at a fixed wavenumber 

- Adding Doubling 

- DISORT 

- Our Way: Parameterization and LUTs, No Adding Doubling 
Calculation, 2-3 orders faster than DISORT. 

• PCRTM - Principal Component-based Radiative Transfer Model 

- Uses PCA to compress information content 

- Reduce RT calculations by at least 3 orders of magnitude 

• Expectation 

- over 4-5 orders faster than a reference radiative transfer model 
such as MODTRAN 
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A Simple Cloud- 
Earth System— 
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M^in Scattering Process: 

Beam beam; beam diffuse; 
diffuse beam; diffuse diffuse 
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(^ Parameters Related to the Main 

Scattering Processes 
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All of these parameters were calculated using 36-stream DISORT to form the LUTs. 
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nsize the Bidirectional Reflection 
Distribution Function 


Rc = R(r,D e ,fj 

• Size of R c LUT in visible range with 100 wavelengths for 10 
different particle size: 

45x10x30x30x30x100x4 = 4.86 GB 


• Size of Q LUT: 

10x30x30x30x100x4 = 108 MB 


R c = R(t, Q{D e , ju 0 ,ju, <j> - A)) 
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ffective Multi-Scattering Stream 
Model (EMSSM) _ 
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fective Multi-Scattering Stream 
Model (EMSSM) 


Radiative Transfer Eqaution: 
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^^Effect ve Multi-Scattering Stream 
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Single Scattering: 

* Source: Directly from Sun, the intensity is known 

* Incident Angle: Solar angle (known) 

Multi-Scattering Stream: 

* Source: diffused sun light by cloud, intensity is unknown, fitting parameter 

* Effective Incident Angle: unknown, fitting parameter 

Each scattering stream may be considered as an effective single scattering with unknown source 
intensity and unknown incident angle. In this work, we using nonlinear regression method to fit the 
BRDF data to get these unknown effective source intensity and unknown effective incident angle 
for each of the scattering streams. These parameters are then used in the fast model calculation. 
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curacy of the Parameterized Ice 
Cloud BRDF @ 600 nm 
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Conditions 

Maximu 
m Error 

Mean 

Error 

Standard 
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Error 

X = 600 nm; 

De = 2.0 - 10.0 pm; 
i = 0 — 550.0; 
SZA = 0 - 90°; 

0.0189 

-4.1 x 
10-e 

1.4 x 10- 3 

VZA = 0 - 60°; 

VAA = 0 - 360° 
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Counts Counts 


r curacy of TOA Reflectance Calculation Using PCRTM Fast 

Model @ 600 nm 
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Results for R s = 0.5 

^gas above — 0.5 . 
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A, = 600 nm; De = 2.0 - 10.0 ^m; x = 0 - 550.0; x gas _ be , ow = 0 - 10.0; SZA = 0 - 90°; VZA = 0 - 60°; 
VAA = 0 - 360° 

All parameters except particle size were selected to be off grid. 
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Counts Counts 



Accuracy of the TOA Reflectance Calculation Using 

PCRTM Fast Model @ 3000 cm -1 
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Results for R s = 
0 . 1 , 

T =01 

u gas_above Vm ' 



u = 3000 cnv 1 ; De = 10.0 - 180.0 pm; x = 0 - 20.0; x gas below = 0 - 5.0; SZA = 0 - 70°; VZA = 0 - 60°; 
VAA = 0 - 360° 

All parameters except particle size were selected to be off grid. 
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BRDF 


Generation of LUTs: Stream Number Dependence in 

DISORT Simulation 
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BRDF BRDF 


ration of LUTs: Remove the Stream-Dependence in DISORT 

Simulation for the LUTs 




0.36% of the data were replaced due to stream-dependent errors. 
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Downsize of LUTs 



• BRDF: 3,300 MB to 23.68 MB 

• T bd : 157 MB to 11.66 MB 

• T DB : 157 MB to 11.66 MB 

• R dd : 5.23 MB 
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BRDF 


Rebuilt of LUTs in PCRTM SOLAR 




CLARREO Science Team Meeting, Lawrence Berkeley National Lab, April 28-30, 2015 ( Xu.Liu-1@nasa.qov ) 


Accuracy of the Rebuilt LUTs in PCRTM_SOLAR 


Accuracy of the Rebuilt BRDF 

LUTs 

(ice cloud) 

(compared to 52-stream DISORT) 
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v = 1800 - 3000 cm 
D = 10 - 180 pm 

x = 0 - 20 @ 640 nm 


SZA = 0 - 85 
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Parameters Used 


Cloud Type 

Ice 

D e : 10-180 pm 

SZA: 0 - 85° 

x cioud : 0 - 20 @ 640 nm 

VZA: 0 - 70° 

v: 1800 - 3000 cm -1 

VAA: 0 - 360° 
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A Reflectance Simulation: PCRTM_SOLAR 

v.s. DISORT 
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In the case of background scattering with VAA = 0° or front scattering 
with VAA = 180°, max errors were obtained. This is due to the 
inaccurate simulation of DISORT at these specific directions 
(singularity issue in DISORT). 
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ccuracy in TOA Reflectance: PCRTM_SOLAR v.s. 

DISORT — 


Difference in Reflectance at TOA 
(PCRTM_SOLAR v.s. 52 Stream 
DISORT) 


|AR| < 10- 3 

93.30% 

10- 3 < |AR| < 10- 2 

6.65% 

|AR| > 10- 2 

0.05% 

Max AR 

0.0337 


Parameters Used 


Cloud Type 

Ice 

D e : 10-180 pm 

SZA: 0 - 85° 

x cloud : 0 - 20 @ 640 
nm 

VZA: 0 - 70° 

v: 1800-3000 cnv 1 

VAA: 0 - 360° 


^below cloud 0-5 


Reasons for the difference: 

1. The stream-dependent errors in DISORT simulation ( - 0.36%). 

2. The inaccuracy in front and back scattering directions in DISORT simulation 
(The percentage with |AR| larger than 10 3 reduced from 6.65% to 5% if these 
cases removed). 

3. The interpolation errors in PCRTM_SOLAR. 
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COMPUTATION SPEED: 

PCRTM SOLAR v.s. PCRTM THERMAL 





Example Satellite Sensor: IASI 0.25 cm 1 Spectral Resolution Full 
Channel Set 


- PCRTM_SOLAR: 4.89 ms/run 

* 1000 runs with the following parameters: 

SZA = 10°, VZA = 60°, VAA = 72.5°, u changes with 439 different 
values, x above changes with wavenumber, x be , ow changes with 
wavenumber, x doud = 1.025, De = 48 jam, Rs = 0.02 

- PCRTM THERMAL: 6.06 ms/run 


* 1000 runs with the following parameters: 

VZA = 60°, u changes with 439 different values, x above changes with 
wavenumber, x ba | 0W changes with wavenumber, x C | 0Ud = 1.025, De = 
48 jam, Rs = 0.02 

- PCRTM SOLAR is a little bit faster than PCRTM THERMAL. 


• Integrate PCRTM_SOLAR to 
PCRTM will NOT influence the 
computation speed of PCRTM 
greatly. 
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COMPUTATION SPEED: 
PCRTM SOLAR v.s. DISORT 



• Example Satellite Sensor: IASI 0.25 cm 1 Spectral Resolution 
Full Channel Set 

- PCRTM_SOLAR: 0.52 ms/run 

• 1000 runs with the following parameters: 

SZA = 10°, VZA = 60°, VAA = 72.5°, u changes with 439 
different values, x above changes with wavenumber, T be , ow constant, 
T C |oud = 1 -025, De = 48 |xm, Rs = 0.02 

- DISORT (52-stream): 63.61 s/run 

• 439 runs with the following parameters: 

SZA = 15°, VZA = 50°, VAA = 30°, u = 2500 cm’ 1 , i above = 0.3, 
T beiow = 4.8, T doud = 1 .5, De = 10 |jm, Rs = 0.2 

• PCRTM_SOLAR is 122,327 times (5 orders) faster than 
DISORT!!! 


CLARREO Science Team Meeting, Lawrence Berkeley National Lab, April 28-30, 2015 ( Xu.Liu-1@nasa.qov ) 


Example for IASI 0.25 cm 1 Spectral Resolution 

Full Channel Set 
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Radiance (mW.m .sr 


Example for IASI 0.25 cm 1 Spectral Resolution 

Full Channel Set 
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Example 


for IASI 0.25 cm 1 Spectral Resolution 
Full Channel Set 
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